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1 Introduction 



In the Standard Model, the couphngs of quarks and H^-boson are given in the form of 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. In this picture, quark mixing 
and CP - violation are closely related to each other. In order to check whether this 
picture really describes CP - violation and, more generally, whether the Standard Model 
gives a proper description of weak decays of hadrons, one must determine precise values 
of the CA'M-matrix elements - \Vij\. One of the most poorly known among them is 
\Vub\- Naturally it is to be extracted from experimental data for the semileptonic b ^ u 
transition. Experimentally, these processes are very difficult to measure, because of the 
dominant b ^ c transition. Thus, b —>■ u are clearly observable in a small fraction of 
the phase space, i.e. beyond the lepton momentum spectrum for b ^ c. Theoretical 
description of the inclusive B Xuii'e decay is tremendously difficult and suffers from 
large uncertainties. On the other side, there is a general hope that \Vub\ can be extracted 
from the exclusive channels, particularly from B — > piui. For reliable extraction of \Vub\, 
we must have precise experimental data on branching ratios as well as a hadronic matrix 
element accurately calculated from QCD. Recently, CLEG collaboration |^ observed 
and measured the branching ratios for this decay. The reported values are: 



R ^rO +P--\-\ (3-88 ± 0.54 ± 1.01)10"^ WSB 
Br{B t Vf)-^ (2.28 ±0.36 ±0.59)10-^ ISGW 



This branching ratio depend upon theoretical input for the efficiency calculation. They 
used the quark model predictions (WSB and ISGW [0). We see that differences 
between two results are due to the models they used, i.e. theoretical uncertainties are 
important source of error. Theoretically, the problem is how to calculate the hadronic 
matrix element and the corresponding weak transition amplitude since it receives hard 
non-perturbative (low-energy) QCD contributions. Unfortunately, up to now there is 
no theoretical tool by which one can calculate the corresponding form factors exactly in 
the whole physical region for this decay. 

Most of the problems with quark model calculations are related to a lack of a fully 
relativistic treatement of quark spins. In fact, they calculate a specific form factor at one 
physical point and then assume the functional dependence on in the whole physical 
region. The most popular assumption is the polelike one. A step forward in calculating 
form factors was the chiral perturbation theory of heavy hadrons, where the chiral and 



heavy quark symmetries were used to construct a phenomenological lagrangian p6|, p7 
However, only the small range of close to q^^x ^e covered, and again the ansatz 
on functional dependence of form factors must be adopted. Different versions of the 
QCD sum rules were employed as well. These calculations allow to conclude about the 
functional dependence, except when close to q^ax- Still, the various QSR give different 
results. Finally, there are lattice QCD simulations which, along with the QCD sum 
rules, are the only methods to treat non-perturbative QCD in a consistent way (the 
nice review of the heavy mesons phenomenology from lattice QCD can be found in 
Ref . pi|] ) . But, due to the t/V^-cutoff, the relevant matrix elements are calculated for 
niQ ~ nic and q^ < m^. HQET is then used to control the extrapolation to nib and 
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consequently the values of accessible are restricted to a small range close to g^ax- 
However, this is the phase-space region which is expected to contribute to a substantial 
fraction of the decay events and is above the endpoint for charm production in the decay. 
For branching ratio, we need the values of the hadronic matrix element at small too. 
Then the corresponding form factors must be extrapolated to this region. To perform 
this, we need an extrapolation law which is always the hypothesis we make. The form 
factors values after extrapolation to small q^ are extrapolation-law dependent. In this 
respect, the situation is not very different from the quark- model approach. Still by 
lattice results we may constrain the unitarity bounds which is the subject of this letter, 
and hence try to reduce the possible choice on the scaling laws for extrapolation to small 

The authors of Ref. ||T6[| proposed another model-independent way to extrapolate \ Vuh\ 
using the double Grinstein type ratio and more specifically the following decay modes: 
B — > plvi, B — > K*!/!?, D —>■ pii>£ and D K*lvi. However, the rare decays B — > K*^!) 
and D — > K*lvi have not been observed and there is much to be learned on form factors. 

In Sec. 2 of this paper, we give the necessary definitions and motivations for this 
analysis. In Sec. 3 we generate the unitarity bounds on the form factors which are 
constrained by lattice results in Sec. 4. In Sec. 5 we examine the functional dependence 
Ai(g^) by additionaly constrained bounds. Concluding remarks are given in Sec. 5 and 
6. 



2 Form factors 



The hadronic matrix element for the B plvi decay is parametrised as: 

2V{q^^ 



<p{p')\J^\B{p)> 



+ 



M + m 

^Jf^{p + pr{< p)-^—i<p)A{q'){p-pr (2) 
M + m q^ 



where = (V^ — A^) = ^7^(1 — 7^)6; M and m are the masses of the B and p mesons, 
p and p' are their momenta respectivelly; q = p — p' is the momentum transferred to the 
leptons and e'" is the polarization vector of the p-meson. The form factor A{q'^) does not 
contribute to the decay rate in the limit of massless leptons. In this case {i.e. i = e, p is 
a very good approximation) , we have: 
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q'X'/\q') [\H^{q')\' + \H_{q')\' + mq')\'] . (3) 



X{t) = {t + — m^)^ — 4M^m^ is the usual triangular function. Hq and H± are the 
helicity amplitudes which come from the longitudinally and transversely polarized p- 
mesons and are given by: 



M + m 



{M^ - m' - q^)Ai{q^) 



{M + m) 



-Mq' 
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H±{q') = iM + m)A,{q')T^^V{q') (4) 

The dynamics of this decay is described by y4i(g^), A2{q'^) and Vlq"^) - Lorentz invariant 
form factors, which are obviously functions of (= t). The physical region < < 
Qmax (20.3 GeV"^) is very large. In dispersion relations approach, the form factor Ai,2(^) 
can be associated to J-^ = 1+, while V{t) to = 1" intermediate state. This suggests 
a nearest pole dominance assumption on the behaviour of the form factors. -ELCf^ and 
APi?0 used this approach to extrapolate to = for all form factors: 

h{<f) = (5) 
~ — 

where = V,Ai,A2. 

In Ref. it was found that Ai{t) decreases with t, while A2{t) moderately in- 
creases. Ref.|0] suggested that V{t) is consistent with the pole dominance, Ai{t) de- 
creases and A2{t) can be fitted with the pole behavior. Light cone sum rules were used 
in Ref.[0 and ||29| and they conclude that V{t) is steeper than pole, Ai{t) increases but 



is flatter than pole and A2{t) is compatible with pole behavior. UKQCD [Q in their 
analysis show that the pole behavior for Ai{t) is preferred (but with nipoie = 7lli GeV). 
Casalbuoni et al. keep the nearest-pole dominance for V{t) and A2{t), while Ai{t) 
is consisted of polar and polynomial terms, so that their effective lagrangian approach 
leads to the so-called 'soft scaling': 



VA2iq' ) ME±^ . AUq^ )^ (6) 

The common choice is the so-called 'hard scaling' which comes from HQET at leading 
order 1/Mh and in infinite heavy quark mass limit: 

V, Milax) ~ ; v4i(g^,J ~ (7) 

These scaling laws (up to 0{M^^) and log corrections) are used by lattice groups for 
the extrapolation from itlq ~ ttlc to mi,. A very nice discussion about the scaling laws 



can be found in Refs. |p7|, p8| . The other way is to extrapolate first to g^ = and then 
to heavy masses. The advantage is evident since the physical region is smaller and the 
results after extrapolation to g^ = are not that much affected by the ansatz assumed 
for the form factors' behaviour. But the HQET scaling laws are valid only for the 
'f — Imax- The way out was recently pointed out in ref.p9| which states that one can 



extract from QCD that all form factors decrease at g^ = with heavy mass as ~ M^^^"^ . 
We hope that this would help us to reduce the errors in future lattice analyses. There 
were also several quark models employed for a prediction oi B ^ pCvi form factors. 
For instance, in the framework of the light-cone formalism, a relativistic treatement of 



spin was proposed in Ref.[15|, but the form factor values were accessible for g < 0. 
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Dispersion formulation of this quark model was used to relate form factors from < 
region to the physical region by performing the analytic continuation [|1^]. A relativistic 
quark model based on the quasipotential approach was discussed in Ref . |ll3| . From this 
small list of results, in spite of the evident progress of QSR and lattice results, we see 
that no definite conclusion on the functional dependence of the form factors can be 
drawn. Let us derive unitarity bounds for V{t) and Ai(t). In this letter, we concentrate 
to these two form factors and more particularly to Ai{t) which dominates the decay 
rate at large t. This sort of analysis was started by authors of Ref. , and applied to 
iriu^ decay in Ref. ||T9[| . In S-physics, Refs. 



18, 20, ra, 24, 25] used it for the heavy- 



K 

to-heavy meson semileptonic decays. The idea to employ the method for heavy-to-light 
decays was first discussed in Ref. [p!7| . So, most of the material discussed in this paper 



can be found in the above references. 



3 Unitarity bounds 

The starting point is the two-point function: 

Ull''^^J d'xe^"^ < 0|T(j;-^(x) Jr'^t(o))|o >= {q^q- - q^g^^)Ti^^{q') + g^^TlY{<f) i' 
In QCD, both sides of this equation satisfy once subtracted dispersion relations: 

„ dU^'iiq^) 1 /-/mn^'^(t) 



For = 0, we are far from the region where the currents can create resonances 
{{rriQ + mq)KQCD « {'mq + niqY + Q^), so that Xtl '^^'^ ^e reliably calculated by 
means of perturbative QCD. The spectral functions Imll^'^ can be obtained from the 
unitarity relation: 

(q^q-'-q^gnirnUl^it + ie) + g''" ImU^^^it + te) = 

\Y.j c^Pr(27r)^%-pr)<0|J(;_^(0)|r><r|J^,^(0)|0> (10) 
r 

where F are all possible hadron states with appropriate quantum numbers, and the 
integration goes over the phase space allowed to each intermediate state. We proceed as 



in Ref.[0. For fx = u, this is the sum of positive definite terms. By concentrating on 
the Bp-intermediate state, we obtain the strict inequalities. From crossing symmetry, 
we know that the Bp vacuum matrix element is described by the same set of the 
form factors as in (2), but in the different region: ((M + mf < q^ < oo), i.e. on the 
cut. Taking the space-space components of the unitarity relation, we obtain the same 
combinations of the form factors as those entering the decay rate expression. So, it 
suffices to take: 

xv,a={xIa{Q') + \-^,x^AQ")) ^ ■ (11) 
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Finally, we have, 

^"•""W^MmW^"'"'!''"*-'^^ '''' 

Imirid) > -^AHt) [2(M + mr\A,(t)p + \H„{t)\'] e{t - (+) (13) 

where t± — {M ± m)^. Inserting these functions in the dispersion relation at — 0: 

Xv,A = - I dt 1^ (14) 







we obtain the set of inequalities: 



1 r.>^w,T...M2 



dt—^\V{t)\' < 1 (15) 



UttxvU Jt+ 

_l_jy^M^nM.?^\H.?}<^ (16) 

For bounds on the form factors, we satisfy ourselves by calculating Xv,a at leading order: 

D\ f\ x'^(l-xf[(D + 2)(mQ^mg)(mQX^mg(l-x)) + 8mgmQ] 
Xv,A = wr, [^--^) / dx- 



4(27r)^/2 ^ (m^x + m2(l-x))4-f 

(17) 

which in our case (m„ = 0) gives: xv — Xa 



4-n-mi 



Actually, we have three inequalities which constrain form factors in the unphysical 
kinematic region. To obtain the bounds on the physically interesting form factors for 
B p^i'i, we perform the conformal mapping: 



1- ,S) 

by which the whole complex i-plane is mapped onto the unit disc \z\ < 1. More 
specifically: t_ < t < t+ is mapped into the segment of the real axis — 1 < 2; < 0, 
while < t < t_ is mapped into < z < 1. In other words, the physical region for 
Bp — > vacuum transition lies on the unit circle, and for B piu^ decay on the right 
segment of the real axis inside the unit circle. In the z - plane, the inequalities (16,17) 
become: 

^/!^l*.(.)/.(.-)P<l. (19) 

Here, we wrote generically fi = V, A\ and their corresponding functions $j; C is the unit 
circle. The procedure for obtaining the functions ^i{z) is well-known and is the solution 
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of Dirichlet's boundary problem 
Solutions are: 



I|, |2^): its value is known on the circle i.e. |$i(e 



ie\\2 



l~^r~ 32M^m^ (1 + 
3nxv (M + m)5(i_^)| 



;Mm 1 + z 



m) 



Z 2 




2y/Mm 1 + z^ 
M + ml- z 

2VMm 1 + -z' 
M + m 1 - 2 



(20) 
(21) 



The functions ^i{z) are analytic everywhere inside the unit circle. From the point of view 
of analyticity, the problems arise with the form factors. All singularities situated above 
threshold t+ can be absorbed in the phase that can be added in redefinition of |$j(e*^)| 
and eventually will not contribute to our bounds. But the singularities inside the gap 
t_ < t < t+ are important. First of all, there are two poles, one at t = (5.32 GeVY 
Zpoiei = —0.1666, contributing to V{t) and the other at t = (5.73 GeVy Zpoie2 = 
—0.3514, contributing to Ai{t). Since we do not know the residua of the form factors 
at these poles, we will simply remove them. This can be achieved by introducing the 
Blaschke factors: 



V 



Zpolel 



1 



zz. 



Pa, 



Zpole2 



polel 



ZZ. 



(22) 



'pole2 



Since the Blaschke factors Py and Pa, are unimodular, after inserting them into (20) 
{^i{z) — > Pi{z)^i{z)) the inequalities remain the same. In fact, this biases our analysis: 
if we knew the values of residua, the form factors would be precisely determined for the 
large values of g^, or at least our bounds would be very narrowed in the whole physical 
region. 

There is also a problem to incorporate subthreshold singularities. They are expected 
in the analysis of the Ai{t) form factor at {Mb** + nm-nY' . So, there are two branch 
points below the threshold {zi = —0.4671,2:2 = —0.7061) . Their effect is negligible in 
our case, as it can be verified by applying the models discussed in Refs.p3 and [ITB . It 
turns out that the bounds on Ai{t) would be relaxed by no more than 1%. The last step 
in deriving the bounds is the construction of the inner product: 



i9i,9j) 



dz 
2'n'iz 



9iiz)9j 



(23) 



and choose gi{z) = ^i{z)Pi{z)fi{z) and g2{z) = (1 — zzl)^^. Then, from the positivity 
of the inner product, we obtain that the determinant of the {gi,gj) matrix is positive, 
i.e.: 



f,{z,)P,{z,)^,iz,) 



f:iz,)^*M)p*iz^) 

1 

l-lziP 



> 0, V^i G IntC 



(24) 
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or 

We see from the picture (Fig.l) that such bounds are not at all restrictive (|V^(t)| < 
16, |v4i(t)| < 18.3). Our nice exercise did not lead to any reasonable restriction on the 
form factors. To constrain them more we can use some of the form factor values that we 
have on our disposal. Similar analysis was performed for the heavy-to-light transition 
but for the case of 5 ^ -kIp^ in Refs.[g2|, Igf. 




Fig.l: Unitarity bounds on the form factors Ai(g^) and V{q'^). 



4 Lattice constrained bounds 

To incorporate n lattice results in our analysis, we define gi{z) = (1 — zz*)~^ 

(i = 1, . . . , n) and construct the matrix (n -|- 2) x (n + 2) whose determinant is positive, 

i.e.: 



1 f:{z)mz)p:{z) f:{z,)^{z,)p:{z,) 

f\{z)P,{z)^,{z) 



f*{z^)<!>;{z^)P:{z^) 



1—zzt 



1—zz* 



> 



fi{Zn)PiiZn)^i{Zn) 



l-z„z* 



(26) 
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For our purpose, we take three (or two) lattice data of the form factors i.e. the 
results for the 5-meson at rest, which are more precise. Since the UKQCD data are 
very accurate, we take their results (all technical details can be found in Ref.]^). In the 
table below, the lattice results of three collaborations are given 0: 



UKQCD i 
24^ X 48 lattice 

/3 = 6.2 
Clover Action 


g2 [GeV^] 


A(g') 

-L \ J / 




V{q^) 

\J- J 


20.3 

15.3t^ 
16.71^ 

14.41^ 


0.461^ 
0.43^2 

0.39l^ 


0.7+? 
0-7-1 


1.6+} 
1.2+} 
1.51} 
1.41^ 


APE 






g2 [GeV^] 




Aiiq') 




v{e) 


18^ X 64 lattice 

(3 = 6.0 
Clover Action 


20.3 
17.6 
16.6 
13.5 


0.43 ± .08 
0.48 ± .16 
0.70 ± .37 
1.06 ± .98 


0.51 ±0.50 
0.55 ±0.70 
1.05 ± 1.10 


1.6 ± .6 
1.2 ±2.2 
1.2 ±4.1 


ELC 


n 


ttice 

[ 

tion 


g2 [GeV^] 


Aiiq') 




V{q^) 


243 X 60 la 
f3 = 6.4 
Wilson Ac 


20.3 
15.8 
12.5 


0.60(6) 
0.53(9) 
0.43(25) 


0.49(22) 
0.29(37) 


0.77(16) 
0.56(20) 



We see from the table, that the lattice results are quite far from being satisfactory when 
a behavior of form factors is to be studied. Results quoted above are concentrated in the 
vicinity of g^a^,, and must be extrapolated to small q^ . With our constrained bounds, we 
want to restrict the values of the form- factors in the region of small and intermediate g^. 
With lattice data incorporated, the constrained bounds are obtained from (27) (thus, 
we take first three/two results from the table above). Resulting bounds are shown on 
the pictures (Fig.2). 

^The Wuppertal group also studied this decay by extrapolating first to = and then to the 
heavy mass. 
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20.0 



Fig. 2: Lattice constrained unitarity hounds on the form factors Ai{q'^) and V{q^) plotted 
with lattice UKQCD results. 



From the pictures, we see immediately how stronger the bounds are (for instance, 
-0.22 < Ai(0) < 1.14 and -0.17 < 1/(0) < 0.72). They are far more restrictive than 
the previous ones but the problem of the form factors behavior remains unsolved. Beside 
statistical, we did not incorporate other errors in this analysis^. This will be done in the 
next section. If we try to see typical ansatze taken for extrapolations to = 0, we see 
that most of them fall inside the region of allowed values for the form factors at hand. 
Still, we can try to modify our analysis in the step (19). 

5 More constrained bounds 

Again, we perform the conformal mapping: 



" V iv(t - 1-) ^'^^ 

where iV is a constant. The endpoints of the physical region are mapped as; 

^ .^..^^-^^^= (28) 



t = Qiax ' ' Z„an = - | | ) (29) 



^The complete treatement of errors is discussed in Ref.p2[ 
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The inequalities that we derived in Sec. 2 remain the same, but the functions (21,22) 
now become: 



= , =- -(l + z)^(l- z)-^ i 1 + ^ I + -77 (30 

V3^(M + m)5^ ^ ^ ^ V ^/N{l + z)J \ M + m 1-zJ ^ ' 

A\/2MmN ,^ , ,3,^ x-5 A , l-z \ ^ , 2^/NMml + z'^ ^ 



Of course, for A?" = 1 we recover (21,22). 

Let us choose in such a way that Zmax = —Zmin- This gives ^ 1.5 and conse- 
quently Zmax = —Zmin = 0.1011. This mcaus that 2; is a small kinematic parameter in 
the whole physical region, and that we can Taylor expand the functions ^i{z)Pi{z)fi{z), 
around z — i.e.: 

^ oo 

and from inequalities (17) we extract the additional constraint: 

oo 

J]|a„|'<l (33) 

n=0 

Actually, the coefficients in the series can only be obtained from the data (In the 
case of S ^ D^*HPe one coefficient is obtained by the help of HQS which gives the 
absolute normalization of the form factors at q'^^ax- case, we do not have such an 

advantage.). It implies that we have to truncate our series. By taking the first k terms, 
and using the Schwartz inequality and condition (34), we can estimate the truncation 
error as: 



Hr\ 



= max \fi{z) - f^{z) I < max ■ 



1 



< max 



P,{z)Mz)\\^^ 
1 

P^izj^MVT^ 



k+l 
k+1 



For the truncation errors in our case, we have: 



(34) 



k 




Atr[Ai{t)] 







1.761 


2.672 


32.226 


1 


0.178 


0.269 


3.267 


2 


0.018 


0.027 


0.330 
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As it was already mentioned, for conservative bounds we take three/two lattice data 
for the form factors Ai{t) [k = 2) / V{t) {k = 1). We display the resulting plot for 
V{t) even though the truncation error is 'large'. In Ref.[^, it was noticed that we can 
incorporate other uncertainties by choosing the constant on the r.h.s. of (34) greater 
than one. Namely, in calculation of Xv,a{0), we may include 0{af) corrections which 
contribute ~ 18% of the one loop contribution. This corrections, and uncertainties on 
nib = 4:.8~^l GeV, as well as contributions of the subthreshold singularities relax the 
constraint (34) by no more than 20%, i.e. Yl'^=o l^nl"^ < 1-2. By this condition, we can 
determine the value of the coefficient 03(02) in the expansion for Ai{z)(y{z)). Since we 
consider the bounds, we take the sign which leads to more conservative bounds. With 
this in mind, we obtain: 





Oo 


Oi 


02 


03 


upper 


0.0178 


-0.0393 


-0.3474 


-1.0379 


lower 


0.0141 


-0.0956 


-0.7279 


0.8129 



V{z) 


Oo 


Oi 


02 


upper 


0.0155 


-0.0412 


-1.0946 


lower 


0.0123 


-0.0589 


1.0938 



Besides statistical errors, in Ref . @] the systematic errors were estimated. Apart from 
quenching, they quote 11% for Ai{t) and 15% for V(t) of systematic errors. We relax 
the resulting bounds by the value of these errors. Final bounds are plotted on Fig.3. 




Fig. 3: More constrained bounds on the form factors ^1(^7^) and V{q'^) (see text) and 
lattice UKQCD results (bounds are relaxed by the value of systematic errors). 



From these figures, we may see first that Ai{t) increases with t. We also notice that 
Ai{0) is small and Ai{0) < 0.18. On Fig.4, we plot some predictions on functional depen- 

^We used also ELC and APE results to generate the bounds. Naturally, these bounds are much 
weaker. 
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dence Ai{t) against our bounds: On Fig. 4c, the comparison with the nearest pole behav- 
lT| shows that Ai must be flatter. Lattice groups {UKQCD, APE and ELC) tried 



lor 



to fit with the pole behavior but with the masses nipoie = 7tl GeV; 5.99±0.62G'e1/; 5.62± 
O.llGeV respectivelly. If we want Ai{t) to fall within our bounds, the pole mass should 
be rupoie > 6.5 GeV which is quite bigger than the nearest pole mass. This can be inter- 
preted as if radial and orbital excitations had more impact on the form factor behavior 



in the physical region. For B iriui in Ref.[30], the authors tried to study such effects 



by taking into account the first few orbital excitations (uncertainties are large since the 
values of fs, are almost unknown). On Fig. 4b, we plot UKQCD prediction by taking 
mpoie = 7 GeV. In Ref . , they obtain that Ai(t) should be the nearest pole plus a 
polynomial term and more specifically, they take a polynomial term to be a constant. 
This is plotted on Fig. 4a and it seems that linear term should be included too. Finally, 
on Fig.4d the functional dependence Ai(t) predicted by light cone sum rules is plotted. 



12 




Fig.4: ajEffective lagrangian (CHPTHH, Ref.j^) prediction of Ai{q'^); h) UKQCD 
pole fit (see text and Ref.j^); a) Nearest pole dominance /|77|/; d) Light-cone QCD sum 
rules prediction 0/ Ai(g^); are plotted against the bounds from Sec. 4 (dotted) and Sec. 5 
(dashed) 



Note that the bounds generated in Sec. 4 are relaxed by the value of the systematic 
errors. 
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6 Conclusion 



In this paper, we studied the unitarity bounds on the form factors for B p£ui decay. 
Form factors and perturbative calculation of a two-point function are related to each 
other by crossing symmetry and a dispersion relation. By conformal mapping we obtain 
the 'unrestrictive' bounds on the physically interesting form factors. The presence of 
poles was taken into account by corresponding Blaschke factors which simply remove 
them, since we do not know the residua i.e. QpBB oicfpoie- The bounds obtained in this way 
are constrained by the lattice results (more specifically, results obtained by UKQCD). 
We took only these values where the heavy meson was at rest, and with errors included 
to make our bounds more conservative. Besides the bounds obtained in this way, we 
wanted to study the functional dependence of the form factors. With more constrained 
analysis, we see that with present data we can study only Ai{t) which dominates the 
decay rate for large values of t. We confirm that Ai{t) cannot decrease, but moderately 
increases, i.e. it is flatter than the nearest pole dominance hypothesis would give. The 
value of Ai(0) bounded in this way, is surprisingly small (^i(O) < 0.18). In these final 
bounds, we did not include the errors of quenching. If we take them into account and 
relaxe our bounds by conservative 10%, we would have Ai(0) < 0.20. For V{t) we can 
not conclude whether it behaves like pole or it is steeper due to the truncation error 
that we have. In fact both, pole and double-pole functional dependence, are fully inside 
the region allowed by our bounds. The analysis performed in this paper is even more 
needed for B ^ K*^ form factors. This work is in progress. 
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